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ABSTRACT 

We show that a blazar classification in BL Lacs and Flat Spectrum Radio Quasars 
may not be adequate when it relies solely on the equivalent widths (EWs) of optical 
lines. In fact, depending on redshift, some strong emission lines can fall in the infrared 
window and be missed. We selected a sample of BL Lacs with firm redshift identifica¬ 
tion and good visibility from Paranal. We targeted with the X-shooter spectrograph 
the five BL Lacs with z > 0.7, i.e., those for which the Her line, one of the strongest 
among blazars, falls outside the optical window and determined the EW of emission 
lines in both the infrared and optical bands. Two out of five sources show an observed 
Ha EW > 5A (one has rest frame EW > 5A) and could be classified as FSRQs by 
one of the classification schemes used in the literature. A third object is border-line 
with an observed EW of 4.4 ± 0.5A. In all these cases Ha is the strongest emission 
line detected. The Ha line of the other two blazars is not detected, but in one case it 
falls in a region strongly contaminated by sky lines and in the other one the spectrum 
is featureless. We conclude that a blazar classification based on EW width only can 
be inaccurate and may lead to an erroneous determination of blazar evolution. This 
effect is more severe for the BL Lac class, since FSRQs can be misclassified as BL Lacs 
especially at high redshifts (z > 0.7), where the latter are extremely rare. 

Key words: BL Lacertae objects: general - quasars: emission lines - radiation mech¬ 
anisms: non-thernral - radio continuum: galaxies 


1 INTRODUCTION 


Blazars are active galactic nuclei (AGN) with strong jets 
forming a small angle w.r.t. the line of sight and emitting 
variable, non-the rmal radiation across the entire electromag¬ 
netic spectrum llBlandford fe Reesl Il978l : lUrrv fe Padovanil 
Il995l) . Although intrinsically rare, blazars are being detected 
in increasingly larger numbers by extragalactic surveys. 
In fact, although on l y 3, 5 00 blazars are currently known 
l|Massaro et al.l .'2Q09. 20 G), th eir number is steadily grow - 
ing thanks to, e.g., the Fermi llAbdo et al 11 201 Ual 20ll]), 
the optical Sloan Digi tal Sky Survey (SPSS: IPlotkin et all 
bold ), and the Planck IPlanck Collaborationll201ll ) surveys. 
Some faint blazars are also b eing detected as serendipitous 
sources in Swift-XRT images liTurrizianill201Cll ). 

Blazars come in two main subclasses, whose major dif¬ 
ference lies in their optical properties: 1) Flat Spectrum Ra¬ 
dio Quasars (FSRQs), which show strong, broad emission 
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lines in their optical spectrum, just like radio quiet QSOs; 
and 2) BL Lacs, which are instead characterized by an op¬ 
tical spectrum, which at most shows weak emission lines, 
sometimes displays absorption features, and in many cases 
can be completely featureless. Historically, the separation 
between BL Lacs and FSRQs has been made at the (rather 
ar bitrary) equivalent w idth (EW) value of 5 A, rest-frame 
bvIStickel et alJ (Il99ll ) and observed-frame bv IStocke et al.l 
dl99ll) . However, the search for a possible bimodal distribu¬ 
tion in the EW of the broad lines of rad i o qua sars has been 
unsuccessful. Indeed, IScarpa fc Falomcd (|l997l ) pointed out 
that radio-selected BL Lacs were, from the point of view of 
the emission line properties, very similar to FSRQs but with 
a stronger continuum. Most BL Lacs selected in the X-ray 
band, on the other hand, had very weak, if any, emission 
lines, and IStocke et all (|l99ll ). when studying the proper¬ 
ties of the X-ray selected Einstein Medium Sensitivity Sur¬ 
vey (EMSS) sample, had to introduce another criterion to 
identify BL Lacs, this time to separate them from normal 
galaxies. This was based on the Ca H&K break, a stellar 
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absorption feature typically found in the spectra of ellipti¬ 
cal galaxies at A ~ 4, 000 A. Given that its value in non¬ 
active ellipticals is ~ 50%. IStocke et ahl d 199 il l chose a max¬ 
imum value of 25% to ensure the presence of a substantial 
non-thermal continuum superposed to t he host galaxy spec¬ 
trum. This was l ater revised to 40% l|Marcha et al.l 1 19961 : 
lLandt et al .1120021 ). 

Based on all the above, it is clear that blazar classifica¬ 
tion depends then on the details of their appearance in the 
optical band where they emit a mix of three types of radi¬ 
ation: 1) a non-thermal, jet-related, component; 2) thermal 
radiation coming from the accretion onto the supermassive 
black hole and from the broad line region (at least in most 
radio selected sources); 3) light from the host (giant ellipti¬ 
cal) galaxy. The strong non-thermal radiation, the only one 
that spans the entire electromagnetic spectrum, is composed 
of two broad humps, a low-energy one attributed to syn¬ 
chrotron radiation, and a high-energy one, usually though t 
to be due to inverse Compton radiation (lAbdo et alj|2010ch . 

Motivated by this observ ational back gr ound , 

Giommi et al.l (201 2. see also iGiommi et al.l 120131 : 

Padovani fe Giommil 1201511 run simulations to try to 
explain the many selection effects that bias the FSRQ/BL 
Lac classification. These simulations included, amongst 
other ingredients, the blazar luminosity function and evolu¬ 
tion, the thermal, non-thermal and host galaxy components, 
and the peak frequency of the synchrotron emission and 
EW distributions. Once the simulated samples had been 
produced, the bias effects were introduced, in order to 
compare the simulations with the real samples. These 
biases are basically the limiting fluxes, to assess whether 
or not a source is detected by a specific survey, and the 
diluted EW, which depends (apart from the intrinsic EW) 
on the shape and intensity of the simulated non-thermal 
emission and the host galaxy light, and determines whether 
a source gets classified as a BL Lac or an FSRQ. One of 
the results of this analysis was that what we call a BL 
Lac can either be an object with a strong non-thermal 
continuum (possibly due to strong Doppler factors), which 
heavily dilutes its lines below the (arbitrary) dividing value 
of EW = 5 A, or a source with intrinsically weak emission 
lines. Radio surveys tend to preferentially select the first 
type of sources (as these are intrinsically radio powerful), 
while X-ray surveys select mainly sources with lower radio 
power, related to low-ionization radio galaxies (LERGs) 
and therefore optically featureless. The consequence of this 
scenario is that the first kind of sources are classified as BL 
Lacs only because they have a strong non-thermal emission, 
or a high Doppler boosting, which swamp their emission 
lines. Such objects are intrinsically FSRQs and so they 
should be classified. On the other hand, most X-ray selected 
BL Lacs are truly (or nearly) featureless a nd constitute the 
genuine BL Lac class IIGiommi e t al .1120121 ). 

The aim of this work is to test these predictions through 
X-shooter optical and infrared spectroscopy of a sample of 
BL Lacs. The paper is organized as follows: in Sect. 2 we 
describe the criteria adopted to build our sample; in Sect. 
3 we deal with the X-shooter observations and the data re¬ 
duction process; in Sect. 4 we present our results; finally in 
Sect. 5 we discuss our findings and draw our conclusions. 


2 SAMPLE SELECTION 


An effective way of searching for FSRQs that have been 
erroneously misclassified is to observe BL Lacs at 2 > 0.7 in 
the infrared band. This is because at these redshifts one of 
the most prominent emission lines in blazars, Ha? (rest-frame 
wavelength A Ha = 6562.8A) moves out of the wavelength 
range where optical spectroscopy is normally done. In this 
scenario, some blazars might have been classified as BL Lacs 
only because one of their strongest and less contaminated 
emission lines falls beyond the spectroscopic range normally 
used to classify blazars. 

To verify this possibility, we selected a small sample 
of BL Lacs to be observed at infrared wavelengt hs as f ol¬ 
l ows: we started from the BZCAT catalogue llMassaro et al.l 
l2009ll and singled out all radio-selected BL Lacs with 
good visibility from Paranal (5 < 25°), z > 0.7, and 
a firm redshift identi ficatio n (at least two lines detected 


111 emission, 


_, Stickel, Fried fe Kiiluj 19931 : iRector fe Stockel 

l200ll : ISbarufatti, Treves fe Falomdl2005lb We then removed 
the objects for which the Ha line was expected to fall in re¬ 
gions dominated by atmospheric absorption features. Since 
the observed redshift distribution of BL Lacs peaks at very 
low redshifts, our final sample includes only five objects, 
which are listed in Table 1. Three of these have a redshift 
z ~ 0.7 — 1, and the Ha line falls in the J band, while the 


other two have a redshift 2 ■ 
H band. 


1.2 — 1.3 and the Ha is in the 


3 X-SHOOTER OBSERVATIONS AND DATA 

REDUCTION 

We observed o ur b l azar sample with X-shooter 
llD’Odorico et al.l 120061 : IVernet et al.l l201lh . a single¬ 
object medium resolution (R = A/AA = 4000-14000) 
echelle spectrograph mounted at the VLT-UT2 telescope. 
The observations were carried out under programme 091.B- 
0092. The choice of X-shooter as operating instrument 
stems from its capabilities to perform simultaneous optical 
and IR spectroscopy with a single exposure. This way we 
can also study the ratio between optical and IR band lines 
without being affected by the variability typical of blazars. 

For both observations the slit width was set to 0.9” in 
the visual (VIS) and near-infrared (NIR) arms, and 1.0” in 
the ultra-violet and blue (UVB) arm. The UVB and VIS 
CCD detectors were rebinned to 1 x 2 pixels (binned in the 
spectral direction but not in the spatial one) to reduce the 
readout noise. With this configuration, the nominal resolu¬ 
tion is different for the three arms: R ~ 5100, 8800, 5300 for 
the UVB, VIS, and NIR arms, respectively. 

All our targets were observed during ESO period P91, 
i.e., between 1 April and 30 September 2013. The list of ob¬ 
servations is reported in Table 1, together with the observing 
dates and the total exposure times. Sources requiring more 
than 1 hour of integration were observed more than once. 
Each observation is constituted by several exposures, taken 
nodding along the slit with an offset of 5" between expo¬ 
sures, following a standard ABBA pattern. 

We processed the spectra using version 1.4.5 of the 
X-shooter data redu ction pipeline ll Goldoni et al.l 120061 : 
iModigliani et af]|201Cf) . The pipeline carries out the follow- 
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iiig steps: the raw frames are first bias-subtracted, and cos¬ 
mic ray hit s are detecte d and removed using the method de¬ 
veloped bv lvan Dokkuml fl2 00il l. The frames are divided by a 
master flat field obtained using day-time flat field exposures 
with halogen lamps. The orders are extracted and rectified in 
wavelength space using a wavelength solution obtained from 
day-time calibration frames. The resulting rectified orders 
are shifted and co-added to obtain the final two-dimensional 
spectrum. In the overlapping regions, orders are merged by 
weighing them using the errors propagated during the re¬ 
duction process. From the resulting two-dimensional merged 
spectrum, a one-dimensional spectrum with the correspond¬ 
ing error file and bad pixel map was extracted at the source 
position. Telluric corrections were applied for spectra whose 
features fall in regions strongly contaminated by sky lines. 
The IRAF0 task telluric was adopted to perform this correc¬ 
tion. Flux calibration is not required by the goal of our work, 
since we are interested only in line spectroscopy. Normalized 
spectra have thus been produced and used throughout this 
work. 


4 RESULTS 
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We computed the rest-frame EWs for the emission lines of 
our sources, which are given in Table 2. The EWs have been 
computed through direct line integration. Errors have been 
estimated resampling the spectra and take into account both 
the normalization uncertainties and the S/N of the spectra. 
We discuss the sources individually in the following sub¬ 
sections. 


4.1 BZBJ0141—0928 (PKS 0138-097) 

For this blazar we detected the Mgll doublet at 2800A and 
three oxygen lines, namely, [Oil] and the [OIII] doublet (see 
Fig. [IJ. A redshift of z = 0.736 ± 0.001 is de rived, in agree¬ 
ment with that reported bv lRector fe Stockd (1200111 . Despite 
a very high S/N ratio, we can not detect Hydrogen emission 
lines, but the Ha feature falls in a region contaminated by 
several sky lines. 

An intervening system is detected at z = 0.5005 ± 
0.0005, featuring absorption by the Mgll AA2796,2803, 
FeIIAA2586,2600 doublets, and the FeIIA2382 line. 


4.2 BZBJ0238+1636 (AO 0235+164) 

This is the blazar featuring the most prominent emission 
lines in our sample. We confir m the r edshift of this source 
to be 0.942 + 0.001 llHealev et al.ll2008h . Five lines have been 
detected, namely, the Mgll doublet, [Oil], the [OIII] doublet 
and Ha. The Hj3 line falls in a telluric region, so only a 
shallow upper limit could be set. The Ha line rest-frame EW 
is higher than 5 A. All the other features detected in the 
optical band have instead values < 5A. BZBJ0238+1636, 
previously classified as a BL Lac due to the EW of its optical 
lines, should be re-classified as an FSRQ according to its 
strong Ha feature in the IR band. The detected emission 
lines are displayed in Fig. [2] 


Figure 1. The emission lines of BZBJ0141—0928. 


iRaiteri et al.l (|2007l 'l have carried out a spectroscopic 
monitoring of this source, finding that the EW of the Mgll 
line was strongly dependent on its magnitude, being smaller 
when the source was brighter. Our rest-frame Mgll EW ~ 
1.5 A is consistent with this trend, as this object was in a 
relatively bright state at the time of our observations (R ~ 


4.3 BZBJ2031 + 1219 (PKS 2029+121) 


This blazar is almost completely featureless. We only re¬ 
port a tentative detection of the Mgll dou blet at the red¬ 
shift z = 1.215, consistent with that by IStickel fe Kiihil 
(111)93 ;. In addition, we detect the inter vening absorber seen 
by Bergeron. Boisse. fe Menardl (1201 if) at z = 1.116, which 
features the MgUA, A2796, 2803 doublet and the [MgIA2852] 
line, which gives more robustness to our redshift determi¬ 
nation. Since the main emission lines should fall in regions 
relatively free of sky lines, this might suggest that we have 
caught BZBJ2031+1218 in a high state, during which the 
thermal emission is overwhelmed by the non-thermal one. If 
this is the case, the emission lines are diluted by the syn¬ 
chrotron emission below the detection threshold of our spec¬ 
trum. Only upper limits for [Oil], the [OIII] doublet, H /3 and 
Ha have been set at the redshift of z = 1.215 (see Table 2). 


4.4 BZBJ2134—0153 (4C-02.81) 

This blazar shows three emission lines, namely, [Oil], 
[OIII]A5007A and Ha. We found a redshift of a = 1.286 ± 
0.001, consistent with that reported by iRector fe Stockel 
(120011) . For the H/3 and [OIIIJA4959A features we could 


1 http://iraf.noao.edu/ 


2 http://www.bu.edu/blazars/VLBAproject.html 
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Table 1. X-shooter observations 


Source 

RA (J2000, hr) 

Dec (J2000, deg) 

redshift 

Observation date 

Exposure time (s) 

S/N 

airmass 

BZBJ0141—0928 

01 41 25.83 

-09 28 42.9 

0 . 733 W 

2013-07-24T08:39:08.945 

2880 

~ 20 

1.09 





2013-08-17T04:57:54.150 

2880 

~ 20 

1.56 





2013-08-19T04:55:01.183 

2880 

~ 20 

1.52 





Total 

8640 

~ 35 


BZBJ0238+1636 

02 38 38.93 

+16 36 59.0 

0 . 94 ( 3 ) 

2013-08-10T09:25:40.669 

1800 

~ 12 

1.34 

BZBJ2031+1219 

20 31 54.99 

+12 19 41.0 

1.215( 2 ) 

2013-08-10T02:38:58.540 

2880 

r^i § 

1.30 





2013-08-10T03:48:49.421 

2880 

r^i § 

1.26 





Total 

5760 

~ 12 


BZBJ2134—0153 

21 34 10.30 

-01 53 17.0 

1.283( 1 ) 

2013-05-27T06:29:07.736 

3600 

~ 15 

1.54 

BZBJ2243—2544 

22 43 26.40 

-25 44 30.0 

0.774( 4 ) 

2013-06-20T08:33:30.168 

1800 

~ 35 

1.02 


(!) iRector fc Stockel teOOlf) : ( 2 ) IStickel fe KmTH dl993h : ( 3 ) iHealev et all (120081 s ) : ( 4 ) IStickel. Fried fe Kuhd (119931 s ) 


Table 2. Rest-frame equivalent widths for our blazar sample. 



BZBJ0141—0928 

BZBJ0238+1636 

BZBJ2031+1219 

BZBJ2134—0153 

BZBJ2243—2544 

Line 

rest-frame EW (A) 

rest-frame EW (A) 

rest-frame EW (A) 

rest-frame EW (A) 

rest-frame EW (A) 

MgIIA2800A 

1.0 ±0.2 

1.5 ±0.4 

Tentative 

< 0.3 

< 0.2 

[OII]A3728A 

0.3 ±0.2 

0.5 ±0.4 

< 0.2 

2.2 ±0.3 

0.6 ±0.2 

H/3 

< 0.2 

< 1.0 

< 0.3 

< 0.4 

< 0.2 

[OIII]A4959A 

0.3 ±0.1 

0.6 ±0.3 

< 0.3 

< 0.3 

0.2 ±0.2 

[OIII]A5007A 

1.0 ±0.2 

1.5 ±0.4 

< 0.3 

2.8 ±0.6 

0.5 ±0.2 

Ha 

< 1.2 

6.7 ±0.4 

< 0.5 

3.2 ±0.4 

2.5 ±0.5 
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set only upper limits. However, the [OIII] doublet falls in a 
region highly contaminated by sky lines. It is interesting to 
note that the H/3 lower limit is about eight times smaller 
than the Ha value, despite the ratio between the two os¬ 
cillator strengths being only ~ 5. We do not detect the 
Mgll emission line reported in iBergeron. Boisse, fe Menardl 
hoill) nor the Mgll intervening absorber at z = 1.2458. 
However, the reported EW of the intervening Mgll system 
(~ 0.2A) is below our detection threshold. The detected 
emission lines are displayed in Fig. [3] 
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Figure 2. The emission lines of BZBJ0238+1636. 


4.5 BZBJ2243—2544 (PKS 2240-260) 

This blazar shows four emission lines, namely [Oil], 
the [OIII] doublet and Ha. The redshift determined 
f or this object is z = 0.780 ± 0.002, consistent with 
IStickel, Fried fe Krilnl (jl993l l. For the H/3 feature only an 
upper limit could be set. Again, we note that this lower 
limit is more than ten times lower than the Ha value, de¬ 
spite the ratio between the two oscillator strengths is just a 
factor of ~ 5. The detected emission lines are displayed in 

Fig- EH 
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Figure 3. The emission lines of BZBJ2134—0153. 
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Figure 4. The emission lines of BZBJ2243—2544. 


5 DISCUSSION AND CONCLUSIONS 

We have looked for Ha emission in the IR spectra of 
five BL Lacs at z > 0.7. We have found that two 

sources (BZBJ0238+1636, rest-frame EW > 5 A and 
BZBJ2134—0153, observed EW > 5 A) could be classified 
as FSRQs by one of the classification criteria that have 
been used in the literature. Indeed, BZBJ0 238+ 1636 was 
defined as an “intruder” BL Lac in iGhisellini et ahl (l201llf . 


This is because its high 7 -ray luminosity and Eddington 
ratio resemble those of typical FSRQs, despite its 7 -ray 
energy index being border-line between the two classes. 
Concerning BZBJ2134—0153, its classification is ambigu¬ 
ous, since i t was consid ered an FSRQ in the 1LAC cat- 
alog (lAbdo et all 2010dll and a BL Lac in the 2FGL one 


llNolan et al. 2012). Its 7 -ray luminosity of ~ 10 47 erg s 1 


and its WISE colors fsee lMassaro et all 1201 ill better place 
this object in the FSRQ class. In addition, given the steep 
NIR-to-UV slope, if the redshift of this object were 2 ~ 2, a 
strong Ly-a line would have been detected, definitely placing 
BZBJ2134—0153 in the FSRQ class. A third source in our 
sample shows border-line (observed) EW (BZBJ2243—2544, 
EW = 4.4 ± 0.5 A). The remaining two objects show very 
weak lines or nearly featureless spectra. 

Our sample shows that, whenever emission lines are 
present, Ha is always the most prominent one (the only 
exception being BZBJ0141—0928, for which it falls in a re¬ 
gion strongly contaminated by sky lines). This appears to 
be the case even if this feature falls relatively close to the 
peak of the synchrotron emission in our sources, which are 
low synchrotron-peaked BL Lacs. I 11 this situation the syn¬ 
chrotron emission, which is non-thermal, should strongly di¬ 
lute the thermal processes such as the Ha emission lines. Ha 
turns then out to be very important if one has to rely on a 
classification scheme that discriminates BL Lacs from FS¬ 
RQs according to the EW of their lines. In this context IR 
spectroscopy turns out to be vital for relatively high redshift 
(2 > 0.7) sources. 

The danger of source misclassification bas ed sole ly on 
EW strength has been recently stressed also bv IRuan et ali 
(|2014h . They studied multi-epoch spectra of more than 300 
blazars, finding that for six of them the EW of their emis¬ 
sion lines crossed the 5 A dividing line in different observa¬ 
tions. These objects have high accretion rates, strong vari¬ 
ability both in the thermal and non-thermal continuum and 
a synchrotr on peak frequency similar to that of FSRQs. 
i Ruan et ali (l2014T l propose that these sources are likely FS¬ 
RQs whose jet axis points extremely close to the line of 
sight. This peculiar geometry causes the emission lines to 
be strongly diluted by the non-thermal continuum when 
the latter is part icularly strong. They conclude, similarly 
to Iciommi et al] (l 2012 l l. that a simple classification based 
on the EW strength may lead to consider these FSRQs as 
BL Lacs and to a wrong evaluation of the BL Lac evolution, 
especially at high redshifts. This effect could be even more 
severe when considering that blazars with redshift larger 
than 0.7, currently classified as BL Lacs, could be in fact 
FSRQs with the Ha emission line in the IR window. Thus, 
a new, EW-independent operational definition for ”BL Lac 
object” is needed. 

Such a picture has also strong implications on blazar 
classification and evolution and on the so-called “blazar se¬ 
quence”, and predicts the association of “real” BL Lacs 
with LERGs and “fake” BL Lacs with high-ionization ra¬ 
dio galaxies (HERGs). In particular, the evolution derived 
for FSRQ samples in the radio, and also 7 -ray, band could be 
biased by the exclusion of relatively higli-redshift misclassi- 
fied sources. Given the role that radio emission is thought to 
play in galax y evolution through the so-called “radio mode” 
accretion ( ICroton et al.|[2006ll , this is also relevant for “AGN 
feedback” via the selection of relatively high-redshift HERGs 
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but especially LERGs, which being intrinsically less powerful 
are harder to identify, already at moderately high redshifts, 
through blazar samples. 
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